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The hydrodynamics of smectic films at an air-water interface is discussed, with particular focus 
on the viscous response of the film under flow normal to the layers. The corrections to the response 
functions of the smectic phase, arising from the coupling between the flow and the smectic order 
parameter, are calculated. The results for the effective viscosity are illustrated by analysing smectic 
film flow in a channel geometry. Two limiting cases of the flow, namely, motion dominated by 
dislocation-induced shear-softening and dislocation-free motion dominated by the permeation mode 
of mass transfer, are studied. The effect of drag from a finite depth liquid subphase is considered. 
£3 ■ The results are compared to those for hexatic and liquid films. 
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I. INTRODUCTION 



The rheological properties of liquid crystalline systems continue to be of considerable interest because of their rich 
and complex behaviour. Thin films of anisotropic molecules, such as Langmuir monolayers at an air-water interface, 
are relevant to many industrial applications, and as such, have been subjected to detailed experimental studies. In 
\ particular, the viscous response of liquid crystalline films is often found to be non-Newtonian. Among the dominant 
£3 causes is coupling of the flow to molecular alignment. In smectic films (crystalline in one dimension, but liquid-like in 
the other), the presence of unbound dislocations becomes a major factor affecting viscous response. When the film is 
riding over another phase, usually water, viscous drag from this subphase, if large, can modify the flow profile of the 
film quite significantly. 

t-H , The coupling between molecular alignment and flow has been seen in some cases in experiments by Mingotaud et 
al. |Q] , Maruyama et al. , and Kurnaz & Schwartz || . The experiments involve Langmuir monolayers of rod-shaped 
molecules that are usually tilted with respect to the surface normal, forming a hexatic phase with anisotropic in-plane 
bond orientations. Typically the film consists of domains of a liquid crystalline phase co-existing with another liquid 
crystalline phase or with the liquid expanded phase (no orientational order). The domains can be distinguished 
through Brewster Angle Microscopy which is sensitive to molecular orientation, making it possible to follow the shape 
and movement of the domains along the flow. There is evidence of nonlinear shear response |^,^| emerging from such 
OS • studies, as well as of the molecular orientation being influenced by flow [|]J|] . 

Shear-thinning has often been observed in experiments involving Langmuir monolayers . A possible explanation 
of this phenomenon is provided by Bruinsma et al. Q in terms of shear-induced defect proliferation. Dislocation 
defects in a solid, if unbound, can relax an applied strain by moving in response to the resulting stress. Since the 
force on such a defect depends on its "charge", oppositely charged defects tend to separate under an external stress. 
Tightly bound pairs cannot contribute to the steady state viscous response, although they can modify the response 
at non-zero frequencies and wave- vectors However, at finite temperatures, bound pairs can dissociate under this 
^ , separating influence which effectively tilts the potential well confining the pair, allowing them to be free beyond a 
potential barrier. Being thermally activated, the free dislocation density, and hence the viscous rsponse, would be 
very sensitive to temperature. Experiments conducted by Schwartz et al. Q on anisotropic hexatic and crystalline 
phases of Langmuir monolayers do indeed see a strong temperature dependence of the critical shear rate for onset of 



?H j non-Newtonian behaviour. 

In this paper we study a simpler problem, the linear hydrodynamics of two-dimensional smectic films in a channel 
flow geometry. Dislocations still play an important role, and it is easier to analyze their effect on the smectic order 
embodied in a single set of Bragg planes. Although we do not study this here, channel flow of two-dimensional 
smectics would also be a promising context in which to explore a tractable model of shear thinning. 

In the absence of external strains, free dislocations or disclinations don't occur in the most ordered two-dimensional 
phases; they are instead bound in pairs of opposite charges by a logarithmic potential. However, in two-dimensional 
layered materials such as smectics or cholcstcrics, there is exponential decay of translational order in both the layering 
direction, and the liquid- like direction along the layers [||. As a result, isolated dislocations have a finite energy and 
exist in a finite concentration at any finite temperature. In these materials, shear response is primarily due to the free 
dislocations; the viscosity diverges inversely as the dislocation density when it becomes small at low temperatures. 
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This divergence is cut off at short length scales by the permeation mode of mass transfer in smectics, where a layer 
distortion induces molecules to jump from layer to layer without affecting the layering structure, allowing the distortion 
to relax over a finite distance. An analogy can be made to the screening effect of supercurrents in a superconductor 
that allow the magnetic field to penetrate a finite distance in from the surface, although the field is expelled from the 
bulk. The "permeation current" plays a similar role with respect to shear in a smectic without dislocations. 

In three dimensions, this analogy has been carried further to predict divergence of the response functions of the 
smectic near the second-order smectic-to-nematic transition. The coupling of the ncmatic order parameter to fluctu- 
ations in the magnitude of the smectic order parameter causes, among other quantities, the permeation constant of 
the smectic and the viscosity denoted 771 in the literature, to diverge However, in two dimensions, a dislocation- 
driven thermodynamic transition to the nematic state occurs at zero temperature dislocations being the analog 
of Abrikosov flux vortices in two-dimensional superconductors. At a finite temperature, the nematic melts into an 
isotropic liquid via a disclination unbinding transition. Below this temperature, local smectic order is disrupted by 
singularities in the phase of the order parameter (i.e., dislocations). However, the local smectic order parameter has 
a finite magnitude, and fluctuations in the magnitude are irrelevant in the renormalization sense JTo| |. Thus renor- 
malization of the elastic coefficients and response functions as in the 3-d case does not occur. The phase of the local 
smectic order parameter is related to the layer displacement, and dislocations, which lead to branch cuts in this phase, 
cost a finite energy as for magnetic vortices in the 3-d superconductor. 

Coupling of the film flow to a subphase (a fluid body supporting the film on its surface) can also significantly 
alter its flow profile. Such experiments have been conducted by Schwartz et al. [ftllPI using Langmuir monolayers on 
water. When the subphase drag dominates the flow, the flow profile becomes semi-elliptical. Stone |l2j has performed 
calculations which confirm this profile and also yield the profiles interpolating between the elliptical and the parabolic 
as the viscosity of the film relative to that of the subphase is increased. The depth of the subphase was also a 
parameter in the calculations, since decreasing the depth results in increased drag. As we show below, Stone's results 
can also be applied to hexatic films. 

The hydrodynamics of partially ordered hexatic films has been studied in detail by Zippelius et al. ]l3| . These 
authors find a correction to the effective viscosity under flow conditions where the hexatic bond orientation is pinned 
at the boundaries, as compared to the case where the bond orientation is free to rotate. The correction comes from 
coupling of the flow to the bond orientation order parameter under the constraint imposed by the boundaries. 

A similar coupling can be enforced by imposing a pressure gradient on the flow. In the experiments conducted by 
Kurnaz & Schwartz || on hexatic film flow, the domain structure of the hexatic mesophase can impose constraints 
on the bond orientation at domain boundaries, thus increasing the viscosity from its bare value. Annealing of the 
domains would then lead to a reduction in the effective viscosity. The experimental signature of this effect would be 
a time-dependent viscous response. Some transients have indeed been observed in these experiments, although other 
factors may be involved, such as domain boundary elasticity |?| and shear thinning. 

To illustrate how to incorporate effects of a subphase into the hydrodynamics of a partially ordered film, we adapt 
the analysis of Ref. |12| to films with hexatic order. In the presence of a subphase, the hydrodynamic equations of 
motion for a hexatic get modified by adding a subphase drag term to the viscous force: denoting by x the co-ordinate 
across the channel, z along the channel, and y along the channel depth (see Fig. Q), we have (with stress matrix Cij), 



d z ir zz + d x a zx + d y a zy 



dg z 
dt 

= 7t' + d x \ vd x g z - ^-d x 8^j - v h d y g z \ fi i m sur face (1.1a) 
86 d x g z 



dt 2p 



+ T 6 K A d^9 (1.1b) 



where ]l3| g is the momentum density, 9 is the hexatic bond orientation order parameter, K A the bond orientation 
stiffness, Tq the corresponding kinetic coefficient, v the dynamic shear viscosity (v = rj/p), Vb the viscosity of the bulk 
subphase, and n' the surface pressure gradient driving the film down the channel. 

Assuming d x 9 is time-independent in the steady state, we have dtd x 9 — 0. Therefore dt9 must be constant acros s 



the channel. Since g z is an even function of x in this flow situation, dt9 is odd, and hence, must be zero. Eq. ( 1.1b ) 
then gives us the coupling between the flow and the bond orientation: 

d 2 x 6 = -—±— d x9z . (1.2) 

2pV 6 K A 



Upon substituting this result into Eq. (1.1a), we find: 
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- ' d x I I v + — — d x g z - Vb d v g z \jn m surface • (1.3) 



dt VV 4 /° r < 

Thus g z obeys an equation of motion identical to that for an isotropic film with a subphase, but with the modified 
viscosity rj + 1/4Tq, and we can take over the results of Ref. As we shall see (Sec. fv|), this simplification does 

not apply to smectic flow. 



In the next section, we briefly review the equilibrium properties of two-dimensional smectic films H. Section [II 



discusses the hydrodynamics of two-dimensional smectics, and the implications of the coupling between the smectic 



order parameter and the flow for the response functions of both quantities. Section IV looks at flow of a smectic film 
in a channel flow geometry, and examines the behaviour in different regimes of the channel width. In Section we 
consider the effect of subphase drag on the smectic film flow as compared to previous results for an isotropic film. 
The results of both sections |^ and [y| are consistent with the results of section |n] for the effective viscosity. The last 
section summarizes the results of this paper. 



II. REVIEW OF SMECTIC FILMS 



Smectics are characterized by a crystal-like periodic modulation of the density along one direction, say, the z- 
direction, and liquid-like correlations perpendicular to it. In two dimensions, we take this to be the x-direction. The 
preferred orientation of the "layers" is also the average direction along which the directors n of the nematic molecules 
are oriented. Although it represents a spontaneously broken rotational symmetry, the layer orientation can be forced 
by boundary conditions on the molecules, or even by flow. Smectic order is characterized by a wavevector qo = z2n/d, 
where d is the layer spacing, usually slightly larger than the molecular length. 

The smectic density wave can be represented as |l4| ] 

p(r)=p (l + V'(r)e lqo - r ). (2.1) 

Here, tp(r) is the complex smectic order parameter: its amplitude represents the strength of the smectic ordering, 
whereas the phase </>(r) = qou(r) describes the phonons associated with broken translational symmetry along the 
layering direction. Phonons in two dimensions are very effective in destroying the one-dimensional translational 
order: The correlation (ip(r)ip* (0)) decays as the exponential of a power of the displacement. Since the square of the 
wavevector appears in the exponent higher harmonics of go m the density modulation are ignored. 
The Landau-Ginzburg free energy takes the form |lq| 



T= d 



^l 2 + \W + |l^| 2 + yip. - iQoSnW + ^-(djn) 2 + ^{d z 5nf 



(2.2) 



Here, K\ and K% are splay and bend elastic constants. The twist elastic constant K2 is absent in two dimensions. 
The coupling between d x ip and Sn = (n — z) • x is required to satisfy the rotational invariance of T . Terms of 
order higher than quadratic in the order parameter and its gradients have been neglected. Well below the mean-field 
smectic-nematic transition temperature, fluctuations in the amplitude of ip can also be ignored, and in the absence 
of singularities in n (disclinations) , Sn can be integrated out. The remaining long wavelength fluctuations can be 
expressed completely in terms of the layer displacement u(r) as ]l4],p| 

T = jd 2 r^B[(d z u) 2 +X 2 (d 2 x u) 2 ], (2.3) 

with B = V>o9o c lh an< ^ ^ 2 = Ki/B- Note that uniform gradients of u along the layer direction (d x u) don't cost any 
energy, because they represent tilting of the layering direction. This important difference compared to two-dimensional 
solids, hexatics, etc. implies that the lowest energy defects in the system, dislocations, have a finite energy Ed and 
are not constrained to be bound in pairs at low temperatures H . 

Whereas a smectic with thermally excited phonons would behave like a nematic with only a splay degree of freedom, 
the presence of dislocations allows for bend in the average layer orientation over scales larger than the typical size £d 
of a correlated "smectic blob" ||, given by £|, = n^ 1 w a 2 D e E,D ^ kBT (ap is a dislocation core diameter, a 2 D ~ dVXd). 
Therefore the long-wavelength behaviour of the smectic is that of a nematic with free energy 



Kx(d x 6N) 2 + K 3 (d z SN) 2 , (2.4) 
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where N denotes the layer normal, and K3 oc As discussed by Nelson & Pelcovits JlB] , non-linearities in the 
nematic free energy modify the nematic Frank constants K\ and K3 such that at scales longer than de^ D ' a the 
nematic can be described by a single Frank constant oc £fj. In practice, this length scale can be very large compared 
to typical system sizes, so one usually sees a 2-Frank constant nematic. A study of th e dynamics of smectic films, 
taking dislocations into account || , yields nematic behaviour corresponding to Eq. J2.4| ) at long length scales, with a 
nematic kinetic coefficient that vanishes like no oc e ~ E o/k B T ^ j ow temperatures. 



III. SMECTIC HYDRODYNAMICS 



The hydrodynamic variables for a two-dimensional smectic are the layer displacement u, and the conserved momen- 
tun densities g x , g z . In this section we focus for simplicity on the dynamics of free-standing smectic films ]l7t , where 
the momentum is conserved to a good approximation. The drag due to a liquid subphase is considered in Sec. We 
assume the smectic to be incompressible, and so neglect density fluctuations, setting the density p = const. 

Consider the stress tensor Uij = Vijkidkgi- From the symmetry properties of the viscosity tensor, it can easily be 
argued that there should be 4 independent viscosity coefficients: r\ xxxx , n zzzz , r\ xxzz — n zzxx , and r\ xzxz — n zxzx — 
Vxzzx = rizxxz- Upon denoting 

h=- 5 -^ = B(d 2 z -X 2 dt)u, (3.1) 
the equations of motion can be written as j[4| 

5 = - + V* (3-2a) 



dt p 

'- - d z p + v z d\g z + vd 2 g z + v ' d z d x g x (3.2b) 



dt 

Jx = - d x p + v x d 2 g x + vd 2 g x + v'd x d z g z (3.2c) 



dt 

where A p is the permeation constant for the smectic, p is the surface pressure, and we have switched to kinematic 
viscosities by dividing by the density, v = n/p, and denoted the four viscosities v x ,v z ,v and v' . The condition of 
constant density: ^| = —d x g x — d z g z = 0, can be used to reduce the number of viscosity coefficients to 3, and decouple 
the (j^-motion from g z . Permeation refers to the dissipative mode of mass transfer in smectics where the molecules 
jump from layer to layer in order to relax a layer distortion. 

Dislocations in the smectic introduce cuts into the displacement field, but it is possible to define the gradient 
s = Vu as a single- valued quantity In the presence of dislocations, Sa; and s z are independent variables, with 
Vxs = — ydm(r), where m(r) is the conserved dislocation density. Now 

| = v| + (i yx JD , (3.3) 

where Jd, the dislocation current, which satisfies 

d t m + V -J D = 0, (3.4) 

is given by 

Jd =n D T i-TT- Vm. (3.5) 



Here we have introduced the kinetic coefficients I 



and the 2-d analog of the Peach-Koehler force 



r, 

L T z 

f = d(Bs z , BX 2 d x s x ); we have also set ks = 1 for convenience. We have imposed the Einstein relation which relates 
the mobility embodied in the first term to the diffusion constant implicit in the second through the common matrix 
r. Since T z and T x correspond to dislocation glide and climb respectively, we expect T z ^> T x . The equations of 
motion are then: 

^ = B(d z s z - X 2 d 3 x s x ) + tt' + vd 2 x g z + v z 8 2 z9z (3.6a) 
-j£ = 8 X ^ + \ p B(d 2 - \ 2 d 4 x )s x + T 2 (n D d 2 B\ 2 dls x - Td z (d x s z - d z s x )) (3.6b) 
-II = d z ^ + X p B(d 2 - \ 2 8 4 x )s z - T x (n D d 2 Bs z - Td x (d x s z - d z s x )) (3.6c) 
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where we have used n' to denote —d z p, and v z 



If conservation of momentum is neglected [j8|, Eqs. (3.6) lead in the limit of long wavelengths and low frequencies 



to a relaxation frequency for s z (which describes layer compression) 

u Sz = -iT x n D d 2 B, (3.7) 
and for s x (i.e., layer undulations) a diffusive frequency 

to Sx (q) = -i ((T z n D d 2 BX 2 )q 2 x + (TT Z + X p B)q 2 ) . (3.8) 



Using the relation Sn — d x u [Q, this last result corresponds to a nematic-like behaviour for the director N = 
z + s x x. Including g z in the hydrodynamic treatment introduces a pair of coupled g-s x modes with both diffusive and 
propagating characteristics. Denoting w 93 (q) = —i(vq 2 + v z q 2 ), the coupled modes have characteristic frequencies 




*(q) « [^f^ ) ± y J + 7^ (3 - 9) 

Propagation dominates for q | x if the dissipation is small enough, which is possible at low temperatures, leading to 

w(q) = ±M\ql (3.10) 

As in the case of hexatics, coupling to the smectic displacement field modifies the viscosity of the film. In the 
absence of dislocations, it is not possible to shear the smectic film without breaking it. The glide motion of dislocations 
facilitates shear deformation. Permeation can also support shear at short length scales. To calculate the effective 
viscosity, we apply an external shear stress to the system, and calculate the steady state response for g z . The 
stress-strain relation calculated in Appendix ^ then yields 



and 



V eff = V ( 1 + v 1 , 2 ) (3-11) 



V f = nJl + — 1— = ) CVI.2! 



VzF x n D d 2 



in the long wavelength limit. This is similar in form to the viscosity correction for hexatics: 77 — )■ 77 ^1 + ^ • 
A similar calculation for s x and s 2 (again using the method sketched in Appendix |a|) shows that 



T z n D d 2 -> T z n D d 2 \1 + — ^ j (3.13) 

and 

Y x n D d 2 T x n D d 2 (l + 1 ) (3.14) 

The permeation constant is not affected by the coupling. 

At low temperatures (or large dislocation energy Ed), no rapidly approaches as e~ E °l T ^ and the effective viscosity 
of the smectic begins to diverge as e ED ^ T , whereas the response of the smectic displacement strains s^ and s 2 to an 
external force goes to the finite value r\~ x instead of vanishing as T z nud 2 . However, since the permeation mode 
relaxes shear over scales shorter than the permeation length S = \J\ p r), the divergence of the shear viscosity is cut 
off for q x » y/?]T z (d/Sl;D) according to 

An £=r 2 > or ' + — 1——Y (3.15) 

T z nDd z + \ v qi V nT z n D d z + 6 z q% J 

Since these hydrodynamic equations are valid only for wavelengths longer than the dislocation correlation length 
in the x-direction, i.e., q x -C ^J 1 where £j_ = (A^) 1 / 3 (|], this rounding off of the viscosity will extend to the 
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hydrodynamic range only if \/i]T z (dA 1 / 3 /6£p 3 ) <C 1. We expect the bare viscosity rj and A = \J K\jB to stay finite 
as T — > 0. However, £d diverges as e Er, ^ 2T . We expect the permeation constant \ p to behave like e~ Ep / T where E p 
is the energy barrier for molecules to jump from one layer to the next. The dislocation kinetic coefficient T z would 
similarly correspond to the activation energy E g for dislocation glide by breaking and reforming of bonds around the 
dislocation core. But this energy barrier should be small compared to that required for molecular hopping across 
the layers, and we shall ignore it in comparison. Then the above condition is satisfied provided -Ed/3 > E p , so that 

£p — > oo faster than S — > 0. 

A similar rounding off is possible in principle for r\ z : 

Vz^Vz + ^— : — 2> ( 3 - 16 ) 

l x n D d z + \ p q% 

however, this saturation is unobservable in the hydrodynamic limit because the dislocation correlation length diverges 
more strongly in the z-direction: q^ 1 3> £|| = (Cd/A) 1 / 3 . Since dislocation climb is similar to the permeation process, 

should behave like X p , and q 2 oc n 4 ^ 3 would be much smaller than ud at low temperatures. 
Although we have assumed the viscosity to be independent of shear rate, at high shear rates we must account 
for shear thinning brought about by the increase in unbound dislocations in the presence of the shear strain. The 
mechanism for dislocation proliferation under a shear stress is similar to that described by Bruinsma et al. j^] for a 
2d crystal of point particles. The stress tilts the effective potential well binding the dislocation pair, allowing the 
pair to dissociate. The extra density of unbound dislocations facilitates further relaxation of the stress so that the 
effective viscosity decreases with increasi ng sh ear ra te (th e shear strain in the steady state depends on the shear rate 



imposed upon the flow). Note from Eqs. (3.11) and (3.12) that the effective viscosities do indeed drop with increasing 
dislocation density nr>. 

The same mechanism would also apply to shear flow in a hexatic film where disclination unbinding would occur in 
the presence of a strain in the bond-orientation angle. Since the orientational order parameter is coupled to the flow 



as in Eqs. (1.1), disclinations can mediate the shear thinning mechanism in the hexatic phase. 



IV. CHANNEL FLOW OF SMECTIC FILMS 



We are interested in flow under shear or a pressure gradient for a film oriented with the layering direction along 
the channel (see Fig. ||). From the previous discussion, we expect a nematic-like profile for s x unless the dislocation 
density is small, in which case we are in the permeation regime and shear is only supported in a boundary layer of 
width 5. We assume the channel is much wider than the dislocation correlation length £d, so that the hydrodynamic 
treatment is valid. Discussion of the effects of a subphase will be deferred to Sec. M. 

In the steady state, we expect s z to be constant, and the equations of motion (|3.6|) reduce to 



% 
dt 

ds x _ 
dt ~ p 



n' + vdlg z - BX 2 8 3 x s x 

X p BX 2 d% + T z n D d 2 BX 2 dls x 



For convenience we reduce the variables by their dimcnsionless counterparts: 



(4.1a) 
(4.1b) 



x — ► ax, g z 



9z, 



B\ 2 



(4.2) 



where 2a is the channel width, and x, g z and s x are now dimcnsionless. We also define the dimensionless dislocation 
density 



A = r)T z riDd 



(4.3) 



and b = a/S, the scaled channel width. In the presence of dislocations, it is convenient to define b' = by/1 + A. 

Upon solving the equations abo ve w ith t he no -slip boundary condition: at x = ±o, g z = and the permeation 
current oc h cx d x s x — (see Eqs. fl3.l| ) and ( 3.2a| )), we find 



1 + A 



(l-z 2 ) 1 f, cosh(6V) 
2 b' 2 \ cosh(6') 



(4.4) 



G 



and s x — d x u where 



u(x) = 



1 



(1-x 2 ) 2 tanhb' (1 - x 2 ) 1 / cosh(6'a;) 
24 + ~lP 2 + \ ~ cosh(6') 



(4.5) 



1 + A 

There are two regimes of interest here: 

• narrow channel : b' <C 1 (i5> ayl + A): 

(l-a; 2 ) 
9z ~* 2 ' 

i.e. we recover the usual Poiseuille profile expected for a structureless fluid. The permeation current X p h has 
the same form, but is smaller than the momentum density by a factor of 0(b 2 ). Also, in the same limit, 

u(x)*^(l-x 2 ) 2 (13-x 2 ). 

Since m«1, the deviations in layer tilt, 9 = —s x , are small. 

• wide channel : b' > 1 (a\/l + A > 6): 

there are two distinct contributions to g z (Fig. ||): 

A \ fl-x 2 \ l - e-^-H) 



1 + Ay V 2 / & 2 (1 + A) 



If A6' 2 ^> 1 (ay/ A(l + A) ^> 5), then the second term can be neglected and dislocations restore a fluid like 
response, but with an effective viscosity 

V eff = V (l + , (4.6) 

confirming the result we found in the previous section. On the other hand, if the dislocation density is so small 
that A6' 2 <C 1 (ay/l + A S> <5 ^ oyA(l + A)), then the second term dominates and the plug flow profile 
characteristic of permeation flow can be seen. 

In the wide channel limit, we also have 

l j A 24 + (1 + A)5/ 2 63 2 + (1 + A)3 fe 4 V 

For the general case, we can estimate the effective viscosity from the flow rate: for Poiseuille flow, the momentum 
flux is given by J^ a g z dx — | - - p . Using this as the definition of rf^ , we find 

rj _ A 3 / tanh6'\ 

rff = TTa + 6' 2 (i + a) [} v~ )' ^ ' 

For b' S> 1 (permeation regime), 

which reduces to ( f4.6| ) for A6' 2 ^> 1, whereas for Ab' 2 <C 1 (and hence A <C 1), 

^ = ^ or V* = ±-, (4.9) 
77 3 3A P 

reminiscent of the result l/A p q 2 we found for low dislocation densities in the previous section. 
For b' <C 1 (6 <C 1) (permeation regime), we have 

which is a small correction due to the permeation boundary layer. 
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V. CHANNEL FLOW WITH A NORMAL FLUID SUBPHASE 



In practice both the film and the barriers that restrict its flow to a channel geometry, float on a volume of 
fluid with viscosity rjb and finite depth H (Fig. Q). If the dimensionless parameter |lj] 

A = r)/rj b H 

is small, then, as for a normal film, the effect of this subphase can be neglected, and the analysis of the previous 
section is sufficient to describe the two-dimensional film flow. If the subphase drag cannot be neglected, the flow 
profile can be calculated through an analysis similar to Ref. [Q. We outline the steps here, and comment on the 
limiting cases. 

We consider a subphase extending from y = (surface with film) to y = —H (bottom). Let v(x,y) be the (z- 
independent) velocity field describing the subphase, v = (v x , v yi v). The velocity profile in the film itself is vo(x) = 
v(x,y = 0). 

In steady state, the equation of motion for v in the bulk of the subphase is: 

(d 2 x +d 2 y )v = 0. 

The boundary conditions on the subphase are: 

• v — for x — > ±oo or y = —H or y = 0, |x| > a, 

• v x and v y — for y — 0,\x\ < a as well. 

We assume that the subphase is incompressible (V • v — 0), which implies that v x and v y must be zero. 
For the film, the equations of motion at the surface are modified by the subphase drag {K\ = BX 2 ): 

tt' + rjdlv - ifi^Sj. - i] b d y v\ y=a = 0, (5.1a) 
d x (v Q - XpKrflsx) + TzUdcPK!^ = 0. (5.1b) 

Once again, we scale variables such that 

x ^ ax, y — > ay, v -> I ^— ] v, and d x s x = ( — ) v p . (5.2) 



as well as 



Ki 



aS, H -> aH. (5.3) 



All quantities are now dimensionless. Since v p is proportional to the "permeation current" , it obeys the same boundary 
conditions as vq. 

The equations of motion can now be written as 

d 2 x v + d 2 v = for - H < y < with v = at y = 0,-H (5.4a) 
1 + (8 2 d 2 - (1 + A)) v p - Adyv\ y=0 = (5.4b) 
9> = (S 2 d 2 x - A)v p (5.4c) 



where 5 and A were defined in Sec. |IV|. Eq. (5.4a) implies that v must have the form 

f°° A(k) 

v(x,y) = / dk — cos kx sinhfc(iJ + y). (5.5) 

J cos kH 

In terms of the Fourier transform A(k), the film velocity 

/>oo 

Vq(x) — / dk A(k)tanh(kH) cos kx. (5-6) 
Jo 

Using Eq. ( 5.4c ), we can express v p in terms of ^o: 

f°° ( k 2 \ 

v p [x) = I dkQ(k)coskx where fi(fe) = A(k)tanh(kH) I -r- I . (5.7) 

7o \ o fc + Ay 
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Upon substituting these relations into Eq. ( 5.4b| ), we obtain a relation for the Fourier transform A{k): 



1 = 



cos /ex. (5.8) 



j dk A(k) (S 2 k 2 + (1 + A)) ( S 2 k 2 + A ) tanh(fcff) + Afc 

The boundary condition vq(x) = for |x| > 1 imposes another constraint on the A(k): 

/ dk A(k)tanh(kH) cos kx = for |x| > 1. (5.9) 
Jo 

This can be satisfied by A(k) of the form 

oo 

A{k)tanh(kH) = k 1 ' 2 ^ 3 a m J 2m _ 1/2+0 {k) (5.10) 



m=0 



where (3 can be chosen for convenience of computation. If this form is substituted into Eq. (5.8), the x-dependence 
can be integrated out to yield an infinite set of linear equations for the coefficients a m : 

a m Gi n (A, 6, A, H) = 2/3 _ x/ "° /2) , n = 0,1,2,... (5.11) 



where 



Gi n (A,S,A,H) = f dfcG(fc;A,,5,A, J ff)fc 1 - 2 ' 3 J 2m _ 1/2+/3 (fc)J 2 „_ 1/2+/3 (fc). (5.12) 
Jo 

The "kernel" for the smectic case, 

G ^ A ' * A ' = fc2 1 1 + aT^) + ^Mury (5 - 13) 

differs from that for a structureless fluid by the term ^Agp (see Fig. ||). This term reflects the correction to rjq 2 we 
found in Section III. As discussed there, the correction is small compared to the normal term k 2 for A> 1, but can 
grow at low temperatures. If A 3> S 2 k 2 , this correction simply appears as an enhancement of the effective viscosity. 
However, when A <C 8 2 k 2 , the correction gives rise to a qualitative change in the velocity profile, characteristic of 
the permeation regime. The plug flow profile (Fig. ^ in this regime is similar to that seen in the case of a thin 
sublayer (H <C 1), with most of the shear occurring in a boundary layer thickness 6 (as opposed to y H/A for the 
"thin sublayer" case). 

When the subphase drag on the film is large (A ^> 1), the profile is very similar to that of a normal film, since 
the second term in the "kernel" , which is independent of the film structure, dominates the flow. The profile is 
semi-elliptical when the subphase is deep [H ^> 1), and resembles plug flow for thin sublayers (H <C 1). 



VI. SUMMARY 



We have studied the hydrodynamics of two-dimensional smectics incorporating dislocations p] in the context of 
shear flow across the layers. The behaviour resembles that of a nematic for length scales beyond the dislocation 
correlation length, with an effective viscosity that represents the role of dislocation motion in making shear possible. 
At smaller length scales, the permeation mode of smectics determines the shear response. These different regimes can 
be observed in channel flow under a pressure head where the channel width sets the observation length scale, provided 
the drag due to the subphase can be neglected. At small dislocation densities, the permeation mode determines the 
flow profile, which evolves from a parabolic profile for channels narrower than the permeation length to a plug-flow 
shape as the channel becomes much wider. In the latter case (the permeation regime), shear is supported only in a 
boundary layer of thickness equal to the permeation length, hence the effective viscosity as determined by the net flow 
rate across the channel grows as the square of the channel width. On the other hand, for large dislocation densities, 
the flow profile is again parabolic, but with the viscosity modified by the dislocation density. 

The dislocation density in turn depends on the shear rate through the shear strain supported by the layers in steady 
state. Under this strain, dislocation pairs in the smectic unbind at a lower energy cost, increasing their equilibrium 
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density and helping to further relax the imposed strain, resulting in a shear thinning effect. This effect has been 
calculated for a 2d crystal of point particles by Bruinsma et al. j5|, and would also be present in the hexatic phase, 
mediated by disclinations rather than dislocations. The flow results in steady state strains in the bond orientation 
order parameter, which are relaxed by disclination motion. It would be interesting to explore this mechanism for 
shear thinning for the smectic films discussed here. 

When the film flows on the surface of a fluid subphase, and drag from the subphase must be taken into account, the 
flow profile depends on the relative viscosities of the film and the subphase as well as on the channel width and the 
subphase depth. The analysis by Stone fll^j , which is supported by experiments, predicts the evolution of the parabolic 
profile into a semi-elliptical or plug-flow profile, depending on whether the drag is due to the subphase viscosity or 
a shallow subphase. In the Introduction, we showed that the same results apply to a hexatic film if described by an 
effective viscosity incorporating the coupling to the bond-orientation order parameter. In the situations described 
above where the subphase drag dominates the flow, these results are also applicable to smectic fillms, since the 
modification to the "flow kernel" of a smectic film with respect to an isotropic film is decoupled from the terms 
describing the influence of the subphase. The subphase drag manifests itself at long length scales where the film 
structure is unimportant. 
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APPENDIX A: CALCULATION OF EFFECTIVE VISCOSITY 



The hydrodynamic equations of motion can be represented schematically as 

- iu;X( q , u) = -L(q)X(q, W ) + f ext (q, w) (Al) 



where X = (<? z ,s x ,s z ), L is a hydrodynamic matrix (see Eq. (|A6|) below), and f ext — (dja^'f, 0, 0), a^f being the 
applied stress. Upon solving for the response to erf J', we find 

X(q, U ) = (-tw + L(q))- 1 /-*(q, to). (A2) 
In the limit uj — > 0, this simplifies to A(q, to = 0) = L(q)^ 1 f ext (q, w — 0), or, 

0,(q,w = 0) = (L(q)- 1 ) f ext (q,uj = 0). (A3) 
Upon inverting this relation, we find 

iq^atf (q,w = 0) = ((W^ g ) ff»(q,w = 0). (A4) 

Upon comparing this result to the definition of the viscosity: <7jj = Vijkidk9h we nn( i the effective viscosity tensor 

v e B{<l). 

Upon writing 



L(qr x 



1 _ DctL(q) 
ai9 J Mmor Sa9i L(q) 

_ T />) L (q)g,s.Minor 3 , Sx L(q) +L(q) gzSz Minor gzSz L(q) 

~ L(q)ff =^ + Mmu,, , L:q: ' (A5) 



we see that the first term simply yields the bare viscosity, and the second is the correction due to the coupling. 

A similar calculation can be carried out for the effcetive response of the displacement gradients s x and s 2 to external 
forces. 

The response matrix, L(q), for a smectic is 
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(vq 2 x + v z ql) -iBX 2 ql 

-iq x /p X p B(q 2 z + X 2 q x ) + T z (n D d 2 BX 2 q 2 x 
-iq z /p -TT x q x q z 



-iBq z 

Tq 2 ) -TT z q z q x 

X p B(q 2 + X 2 q 4 x ) + T x (n D d 2 B + Tq 2 x ) 



(A6) 



We take q z = when calculating , and q x = in calculating , which leads to Eqs. ( 3.11 ) and ( 3.12 ). 
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FIG. 1. Channel flow geometry with a subphase. The film and the barriers forming the channel are on the surface of a water 
table. A surface pressure gradient is applied to the film along the channel. On the right are the typical flow profiles at various 
depths. 



FIG. 2. Channel flow under the influence of a surface pressure gradient along the channel: the smectic layers are normal to 
the flow. 



FIG. 3. In the wide channel/macroscopic limit, there are two contributions to flow: solid line represents the parabolic profile 
due to dislocation-assisted shear, whereas dotted line represents the plug-flow profile characteristic of flow in the permeation 
regime. 



FIG. 4. G s 



-a/G 



isotropic 



as a function of the wave-number k (solid line). G(k) (see Eq. ( 5.12 ) 



is the "kernel" that 
The channel width 



determines the flow profile of the film on a fluid subphase in the channel flow geometry shown in Fig, 
is scaled to 1. For large k, the structure of the film is important, and the ratio differs from 1 by the Lorentzian 1/(A + 5 2 k 2 j 
(dotted line). For small k, the subphase drag (characterized by A and H) dominates the flow, and the ratio approaches 1, 
leading to the hump at small k. 
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FIG. 2. 



velocity profile 




FIG. 3. 




FIG. 4. 



